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ABSTRACT 

Context. The 30 Dor region in the Large Magellanic Cloud (LMC) is the most vigorous star-forming region in the Local Group. 

OStar formation in this region is taking place in low-metallicity molecular gas that is exposed to an extreme far-ultraviolet (FUV) 
radiation field powered by the massive compact star cluster R136. 30 Dor is therefore ideally suited to study the conditions in which 

■ stars formed at earlier cosmological times. 

Aims. Observations of (sub)mm and far-infrared (FIR) spectral lines of the main carbon-carrying species, CO, [C i] and [C n], which 
I | originate in the surface layers of molecular clouds illuminated by the FUV radiation of young stars, can be used to constrain the 

O . physical and chemical state of the star-forming ISM. 

■ Methods. We used the NANTEN2 telescope to obtain high-angular resolution observations of the 12 CO 7 = 4 — > 3, 7 = 7— > 6, and 
13 CO 7 = 4 — > 3 rotational lines and [Ci] 3 Pi- 3 Po and 3 P->- 3 Pi fine-structure submillimeter transitions in 30Dor-10, the brightest 
CO and FIR-emitting cloud at the center of the 30 Dor region. We derived the physical and chemical properties of the low-metallicity 
molecular gas using an excitation/radiative transfer code and found a self-consistent solution of the chemistry and thermal balance of 
the gas in the framework of a clumpy cloud PDR model. We compared the derived properties with those in the N159W region, which 
is exposed to a more moderate far-ultraviolet radiation field compared with 30Dor-10, but has similar metallicity. We also combined 
our CO detections with previously observed low— 7 CO transitions to derive the CO spectral-line energy distribution in 30 Dor-10 and 
N159W. 

Results. The separate excitation analysis of the submm CO lines and the neutral carbon fine structure lines shows that the mid-7 
CO and [C i]-emitting gas in the 30 Dor- 1 region has a temperature of about 1 60 K and a IT density of about 1 4 cirT 3 . We find that 
the molecular gas in 30 Dor-10 is warmer and has a lower beam filling factor compared to that of N159W, which might be a result 
, of the effect of a strong FUV radiation field heating and disrupting the low-metallicity molecular gas. We use a clumpy PDR model 

(including the [C n] line intensity reported in the literature) to constrain the FUV intensity to about xa ~ 3100 and an average total H 
£N| ' density of the clump ensemble of about 10 5 citT 3 in 30 Dor-10. 
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5_j , 1. Introduction low-metallicity gas (increased [Cn]/CO and [Ci]/CO intensity 

■ ratios) are a result of the combined effect of a strong FUV ra- 

The radiative feedback of star formation on its progenitor molec- diation field and of reduce d metallicity. It is therefore difficult 

ular gas is of great interest in modern astrophysics because it to study sep arately the influence of these two parameters on the 

can be a determining factor in the regulation of star formation in p roper ties of the low-metallicity molecular gas. 
galaxies. In particular, there is special interest in the interaction 

of stellar far-ultraviolet (FUV) photons with a lower metallic- The Lar 8 e Magellanic Cloud is an ideal extragalactic labora- 

ity, lower dust-to-gas ratio gas, because it is thought to be the tor y in which t0 stud y star fo rmation in lo w-metallicity environ- 

material from which stars formed at earlier cosmological times. ments - Its proxi mity (50kpc;|_F east 1999|and nearly face-on ori- 



The region of the interstellar medium where the chemistry entatlon (i=35°;| van derMarel & Cioni|200 lb allow us to isolate 

and thermal balance is dominated by the influence of stellar FUV individual clouds and study them at the highest spatial resolu- 

photons is called a p hoton-dominated region (or photodissocia- tlon P osslble for external 8 alaxies - The LMC also shows stmn ^ 

tion region, or PDR; iHollenbach & Tielens 1999 and references var y in 8 FUV radiation fields and therefore very different phys- 

therein). Owing to the reduced attenuation of FUV photons by icalconditions in a low-metallicity (Z Q - 0.4 Z ; | Westerlund | 

dust in low-metallicity PDRs, the CO molecule is more readily UJ29J) interstellar medium. 

photo-dissociated, increasing the abundance of C + and C rel- Regions of special interest are the molecular clouds as- 

ative to that of CO. Therefore, the observational signatures of sociated with the 30 Doradus Hn complex. This complex is 

powered by one of the most dense and massive star clusters 

Send offprint requests to: J.L.Pineda e-mail: in the local universe, R136. Therefore, the 30 Doradus re- 
Jorge. Pineda@jpl.nasa.gov gion exemplifies the special case in which the low-metallicity 
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molecular gas is exposed to a very strong FUV radiation field. 
Among se veral CO-clouds locat ed around R136, the 30Dor- 
10 cloud dJohansson et al.|[l998l) is associated with the center 
of the 30 Dor nebula. The Spitzer 8/mi and CO maps shown in 
Figure Q] show that this region is the brightest CO cloud at the 
center of 30 Dor and is located at a projected distance of 20 pc 
from the center of R136. This proximity to R136 and the pres- 
ence of associated 8/mi sources, which might indicate embed- 
ded star formation, suggest that the gas in 30Dor-10 is exposed 
to an extreme FUV environme nt. The CO peak of 3 0Dor-10 
also associated wit h the [Cn] dPoglitsch et al.lll995l) and FIR 
dWerneret al.ll 19781) peaks. Although 30Dor-10 is the CO peak 
at the center of the 30 Dor nebula, the sub-mm CO and [Ci] 
line intensities are weak er compared with o ther active regions in 
the LMC (e.g. N159W: Ipjneda et al]|2008h . On scales of about 
lOpc, several low-/ t ransitions of 12 CO and 13 CO have been 
observed in this cloud dJohansson et al.; 1998; Mi namidani et al.l 
2011b as well as the [Cn] 158//m line dBoreiko & Betzl 
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19911: iPoglitsch etalJU995h . At larger spatial scales of about 
50 pc, the 12 CO / = 7 -> 6, 12 CO J = 4 -» 3 and [Ci] 
3 Pi — > 3 Pp lines have been ob served with the AST/RO telescope 
(IStark et al.lflWfllKrmll2006l) . 

In this paper we present observations of the 12 CO / = 4 — » 3, 
J = 7 -> 6, and 13 CO / = 4 -> 3 rotational and [Ci] 3 Pj -^ 3 P 
and 3 P2 — » 3 Pi fine-structure transitions in the 30Dor-10 cloud 
in the LMC using the NANTEN2 telescope. The angular resolu- 
tion of these observations allow us for the first time to study sub- 
millimeter lines at spatial scales of about 10 pc. We complement 
this data set with low-/ 12 CO and 13 CO transitions observed 
with Mopra, and ASTE, with similar angular resolution. 

The physical conditions of the molecular gas in 30Dor-10 
are compared with those derived in a region with a similar metal- 
licity but with a more moderate incident F UV radiation field, the 
LMC N159W region dPineda et al.ll2008l) . N159W is located in 
one the edges of the N159 Hn region and is actively forming 
massive stars, as suggested by the presence of seve ral embed- 
ded massive young stellar objects (IChen et al 1 120101) . However, 
the FUV and Ha fluxes in N159 are 10-20 and 7 times lower 
than in the 30 Dor region, respectively, suggesting that although 
N159W is an active region it does not have the extreme FUV 
environment to which 30Dor-10 is exposed. The comparison 
between the physical conditions of the gas in 30Dor-10 and 
N159W will help us to understand the effects of the strength of 
the FUV field on the properties of the low-metallicity molecular 
gas in the LMC. 

This paper is organized as follows. We describe the 
NANTEN2 observations in Section [2] and the observational re- 
sults in Section [3] We derive physical properties of the low- 
metallicity molecular gas using a radiative transfer code inde- 
pendently for the CO and [C i] lines in Section|4] We give a self- 
consistent solution of the chemistry and thermal balance of the 
gas using a clumpy photon-dominated region (PDR) model in 
Section We compare our results to those constrained in the 
LMC-N159W region in Section [6TT] and we study the CO spec- 
tral line distribution in both regions in Section 16.21 We summa- 
rize the results in Section [7] 

2. Observations 

We used the new NANTEN2 4-m telescope situated at 4865 m 
altitude at Pampa la Bola in northern Chile to observe the 12 CO 
j = 4 -> 3 (461.0408 GHz), 7 = 7^ 6 (806.6517 GHz), 
and 13 CO / = 4 -> 3 (440.7654GHz) rotational and [Ci] 
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3 Po (492.1607GHz) and 3 P 2 - 3 Pi (809.3446GHz) fine- 



structure transitions toward LMC 30Dor-10. The line parame- 
ters derived from the NANTEN2 observations of LMC 30Dor- 
10 are listed in Table Q] and the observed spectra are shown in 
Figure Q] 

The observations were made towar d the peak 12 CO J = 
1 — > intensity position in 30Dor-10 dJohansson et al.|[l998l) . 
located at a = 5 h 38 m 48 s .56 and 5 = -69°04'43'.'2 (J2000), us- 
ing the total-power observing mode. We used a reference posi- 
tion at a = 5 h 39 m 33 s 7 and 6 = -69°04'0070 (J2000), which is 
free of 12 CO / = 1 — > emission. The duration of each beam 
switch cycle was varied between 20 and 30 seconds, depending 
on the atmospheric stability during the observations. The final 
spectra result from total on-position integration times of 40 min 
to about 3 hours. The pointing was checked regularly on Jupiter, 
IRC+10216, and IRc2 in OrionA. The applied corrections were 
always smaller than 10". We do not expect significant uncertain- 
ties in the line fluxes due to pointing errors because the CO flux 
distr ibution in 30Dor-10 is smooth at scales smaller than 10" 
(e.g. lMinamidani et al.ll20 08S!). The angular resolution of the ob- 
servations is 38" for the 460^90GHz and 26" for the 810 GHz 
channel, respectively. The main beam efficiencies are 0.5 and 0.4 
for the 460-490 GHz and 8 10 GHz channels, respectively. 

The observations were conducted with a dual-channel 
460/8 10 GHz receiver installed for verifying the submillimiter 
performance of the telescope. Double-sideband (DSB) receiver 
temperatures were ~250 K in the lower channel and ~ 750 K in 
the upper channel. The intermediate frequencies (IF) are 4 GHz 
and 1.5 GHz, respectively. The latter IF allows simultaneous ob- 
servations of the 12 CO / = 7 — > 6 line in the lower and of the 
[Ci] 3 P2- 3 Pi line in the upper sideband. These two lines were 
observed simultaneously with one of the lines in the 460 GHz 
channel. As backends we used two acousto-optical spectrome- 
ters (AOS) with a bandwidth of 1 GHz and a channel resolution 
of 0.37 km s" 1 at 4 60 GHz and . 21km s 1 at 806 GHz. 

As discussed in lPineda et alj {2008), the main sources of un- 
certainty in the absolute line calibration are the precision with 
which the beam efficiencies are determined and the accuracy of 
the atmospheric calibration. We estimate that both sources of un- 
certainty contribute to a total absolute calibration uncertainty of 
20%. In our analysis we added these errors quadratically with 
the radiometric noise, assumed to be given by the formal er- 
rors of the Gaussian fits. Propagating these absolute calibration 
uncertainties for line ratios between lines in the 460-490 and 
8 10 GHz channels, we obtain a relative calibration uncertainty 
of 28%. Line ratios between lines located in the same frequency 
channel do not suffer from the uncertainty on the determination 
of the beam efficiencies (~10%) and we propagate the errors 
only considering the uncertainties on the atmospheric calibra- 
tion (-17%). 



3. Observational results 

3.1. NANTEN2 data 

The spectra observed toward 30Dor-10 are shown in Fig. [T] to- 
gether with Gaussian fits whose results are listed in Table Q] We 
determined the line center and width from the high signal-to- 
noise 12 CO / = 4 — » 3 spectrum and we held them fixed for the 
Gaussian fit of the line amplitude of the other lines. This method 
is adequate because their profiles are, within their limited signal- 
to-noise ratio, consistent with that of the 12 CO J = 4 — > 3 line. 

To calculate intensity ratios between lines in the upper and 
lower frequency channels we need to account for the differ- 
ent beam sizes (38"and 26", respectively). We therefore de- 
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LMC 30Dor-10 NANTEN2/KOSMA 




5 h 39 m 30 3 39 m 00 s 38 m 30 s 

Right Ascension (J2000) 

Fig. 1. {Left panel) Spitze r 8jum emission m ap of the 30 Dor region (Meixner et al. 2006) overlaid by an integrated intensity map of 
the 12 CO ./ = 1 -» line dWonget al.ll201ll) . The contour levels correspond to 20% to 80% of the peak intensity (1 1.3 K km/sec), 
in steps of 10%. The angular resolution of the Spitzer 8/vm and CO maps are 2" and 45", respectively. The center of the R136 
super star cluster and the NANTEN2 38" beam are marked by a cross and a circle, respectively. (Right panel) Observed lines and 
corresponding Gaussian fits. 

Table 1. Gaussian-fit line parameters derived from the NANTEN2 observations of LMC 30Dor-10 and from the literature 



Line 


Amplitude'"' 


Center 16 ' 


FWHM (f " 


Integrated Intensity 




[K] 


[kms" 1 ] 


[kms" 1 ] 


[Kkms- 1 ] 


NANTEN2 data 










12 CO J = 4 -> 3 


5.06 (0.19) 


250.1 (0.053) 


9.02 (0.13) 


48.5(1.9) 


12 CO J = 1 -> 6 


2.89 (0.19) 






27.7 (1.7) 


13 CO J = 4 -> 3 


0.76 (0.14) 






7.3 (1.3) 


[Cl] 3 P! ^ 3 P 


0.58 (0.14) 






5.6(1.3) 


[Ci] 3 P 2 -^P, 


0.58(0.19) 






5.6(1.8) 


Literature data' 1 '' 










12 CO J = 1 -» 


1.1 (0.03) 


250.3 (0.1) 


8.79 (0.23) 


11.3 (0.27) 


12 CO J = 3 -» 2 


5.2 (0.4) 


250.4 (0.1) 


8.4 (0.23) 


49.6 (0.27) 


13 CO J = 1 -> 


0.17(0.02) 


250.7 (0.27) 


5.35 (0.64) 


0.88 (0.12) 


[Cll] 2 P 3/2 ^ 2 P,/2 


6.3 (0.6) 


247.8 (1.6) 


20.5 (1.8) 


109.4 (8.5) 



<n) the rms noise for lines observed with NANTEN2 are for channel widths of 0.37 km s~' and 0.21 km s~ 
for the 460-490 and 810 GHz channels, respectively. 

(6) center and width fixed to 12 CO J = 4 — > 3 best-fit values for NANTEN2 data other than this line 



references for literature data are presented in Section 3.2 
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termined the source extent us ing lower resolution data. As de- the source size from the observed brightness ratio in two beam 
scribed in lPineda et al.l (120081) . assuming a Gaussian source dis- sizes: 
tribution with FWHM S and a source peak brightness tem- 
perature of Ts^eak, the beam filling correction gives a main 
beam brightness temperature in a beam of FWHM &i, of T m t = 
0j/(®j + ©fyTsjeak . The ratio of the main beam brightness 
temperature in two different beams 1 and 2 is given by R \ 2 = 
T mbA /T mba = (® 2 b2 + ® 2 <)/(&l , + With this, we can derive R 2& 2 _ & 2 

®l= • ' ' ■ (1) 

I - "1,2 
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Table 2. Observed line intensity ratios, corrected for beam coupling to an assumed common FWHM source size of 90 





Species 


Ratio 1 "' 


fit error 


calibration error 


total error 


I2 CO J = 


7^6/ 12 C07 = 4^3 


0.52 


0.04 


28% 


0.15 


[Ci] 3 P 2 


-^ 3 Pi / [Cl] 3 P! ^ 3 P 


0.92 


0.37 


28% 


0.45 


13 CO J = 


4^3/ 12 CO 7 = 4^3 


0.15 


0.03 


17% 


0.05 


[Ci] 3 P 2 


^ 3 P, / CO J = 7 -> 6 


0.20 


0.07 


17 


0.09 


12 coy = 


1 -^o/ 12 coy = 4->3 


0.22 


0.05 


28% 


0.08 


13 CO J = 


i -^o/ 12 co/ = i ->o 


0.08 


0.01 


28% 


0.02 


[Cn] 2 P 3/2 


^ 2 P 1/2 / 12 CO J = 4 -> 3 


1.24 


0.13 


28% 


0.37 



ratio of amplitudes from Gaussian fits 



With the source size we can correct the observed intensity ratios 
for the beam coupling of the different beam sizes at two observ- 
ing bands x and y to the source intrinsic intensity ratio R s via 



- mb,. 



T mb , y + 



(2) 



We used observations of the [Ci] 3 Pi — > 3 Po line made with the 
SWAS satellite (Frank Bensch, priv. comm.) with an angular 
resolution of 258". The integrated line intensity of this line ob- 
served by SWAS is 0.44 Kkms" 1 . With NANTEN2 we obtained 
an integrated line intensity of 5.58 K km s _1 in a 38" beam. Thus, 
from Equation (Q]i we obtain a source size of ® s =64", which ac- 
cording to Equation ©, corresponds to a 18% correction to the 
line ratios between the upper and lower channels. We also deter- 
mined the source size by considering the 12 CO J = 1 — » in- 
tegrated line inten sity observed by the NANTEN telescope (e.g. 
iFukui et al.l [20081) of 4.1 K kms 1 in a 158" beam and by the 
Mopra telescope (Jiirgen Ott, priv. comm.) of 11.3 K kms 1 in 
a 33" beam. Equation ([TJ results in a source size of © 5 =110", 
corresponding to a 8% correction. We assumed an intermediate 
value for the source size of ® S -9Q" corresponding to a correc- 
tion of 10%. A source size of 90" is consistent with the angular 
size of 30Dor-10's CO-emitting region shown in Figure [TJ We 
note that assuming a source size introduces an uncertainty of 
about 5% to the line ratios between the upper and lower chan- 
nels. The correction also assumes that the source structure is 
the same at the spatial scales observed with the 460 and 810 
GHz channels. The validity of this assumption likely dominates 
the uncertainties in the line ratios. However, as we will see in 
Sections |4] and |5] the good fit obtained for the physical parame- 
ters supports this assumption. Table |2]shows the line ratios thus 
corrected, which we used for the excitation analysis discussed 

below. 

iKim et alj d2005h presented maps of the 12 CO / = 4 -> 3 
line in the 30 Dor region with a 109" beam using the AST/RO 
telescope. In addition to the 12 CO / = 4 — > 3 component at 
~250kms _1 they detected another component at ~302kms~ 1 , a 
factor of ~3 weaker. We did not find this additional component 
in our 12 CO J — 4 — > 3 spectrum. Considering their derived 
FWHM line-width for this component of 5.1 kms , we derive 
an 3cr upper limit of 0.52 Kkms -1 , which is a factor of ~10 
weaker than the main component of 12 CO J — 4 — » 3 in 30 Dor- 
10. 



3.2. Additional lines from the literature 

To better constrain the physical parameters of the emitting gas 
in 30Dor-10, we combined the lines observed with NANTEN2 
with additional observations of the 12 CO 7=1— >0,/ = 3— >2, 



13 CO J = 1 -> and [Cn] 158 /mi lines. The 12 CO and 13 CO 
/ = 1 — > lines were observed with the ATNF Mopra 22 m tele- 
scope (Jiirgen Ott, priv. comm.) at the position of the NANTEN2 
observations. The data have an angular resolution of 33" and 
a channel width of 0.7kms~'. The 12 CO J — 3 — > 2 obser- 
vations w ere observed with t he AS TE telescope and were pre- 
sented by Min amidani et al.l ( 120081) . The angular resolution for 
both lines is about 22"and the channel width is 0.45 km s _1 . 
The Kuiper Airborne Observatory (KAO) observed velocity- 
resolved [Cn] 158 yum emission toward 30Do r-10 with an an- 
gular resolution of 43" dBoreiko & Be7zlll991l) . The pointing of 
the KAO observations is about 20" farther southwest than the 
NANTEN2 position. We adopted their values without further 
corrections. Indeed, the integrated intensity of the [C n] line from 
the KAO observations of 7.7 x 10~ 4 erg crrT 2 s _1 sr~' is con- 
sistent to within 30% with the integrated intensity at the posi- 
tio n of the NANTEN2 o bservations in the [C n] map published 
bv lPoglitsch et al.ld 19951) . observed with low-velocity resolution 
(Av = 75 km s _1 ) using the FIFI spectral imaging instrument 
on the KAO. In Table [3] we present a summary of rest frequen- 
cies, energies above the ground level, and critical densities of the 
spectral lines used in our analysis. 

We assumed that the uncertainty of the absolute calibration 
of all independently observed lines is about 20%, and propagated 
these errors to obtain a line ratio calibration accuracy of 28%. 
Table |2]quotes the derived total errors, obtained by quadratically 
summing the formal fit error and the calibration uncertainty. 



4. Excitation analysis 

We first constrained the properties of the emitting gas by com- 
paring line intensity ratios and absolute intensities with the re- 
sults of a radiative transfer code. We analyzed [C i] and CO sep- 
arately, but assume that they originate from the same regions 
with identical beam filling factors. This assumption is justified 
considering that the CO and [C i] emission arise from the FUV 
illuminated surfaces of clouds that are unresolved at the scale of 
our observations. 

We used the RADE X non-LTE radiative transfer code by 
Ivan der Tak et al. I (120071) . using the uniform sphere approxima- 
tion, to calculate line intensities as a function of the kinetic tem- 
perature, H2 volume density, and the column density per veloc- 
ity interval N/Av of the species of interest. The collision cross 
sections were taken from the Leiden Atomic and Molecular 
Database (LAMDA; ISchoier et alll2005l) . The code provides the 
average line intensity of a spherical cloud. We therefore scaled 
the data by a filling factor to compare them with the observed 
absolute intensities. As mentioned above, considering a medium 
consisting of clumps that are much smaller than the scale of our 
observations, the C 1 and CO species can be considered as oc- 
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Fig. 2. Comparison between the observations and the predictions 
from the radiative transfer model for the [C i] 3 P2 — » 3 Pi/[Ci] 
3 Pi -> 3 P line ratio (black-solid lines) and the [Ci] 3 Pi -> 3 P 
absolute intensity (gray-dashed lines) for kinetic temperatures 
of 80, 120, and 160 K. The observed [Ci] 3 P 2 -> 3 Pi/[Ci] 
3 Pi — » 3 P() line ratio of 0.92 is shown as a red line with cir- 
cles. The observed [C i] 3 Pi — > 3 Po line intensity of 0.58 K and its 
value scaled by a beam filling factor of 1/8 and 1/16 are shown 
as blue lines with boxes. 




log [n(H 2 )/cm -J ] 



Fig. 3. Comparison between the observations and the predictions 
from the radiative transfer model for the 12 CO 7 — 7 — > 6 to 
7 = 4 -> 3 (black-solid contours) and 13 CO to 12 CO 7 = 4^ 
3 (gray-dashed contours) line ratios for kinetic temperatures of 
80, 120, and 160 K. The plots assume a fractional abundance of 
12 CO/ 13 CO of 35. The observed 12 CO 7 = 7^6 to 7 = 4^3 
ratio of 0.52 is shown as a blue line with boxes while the 13 CO to 
12 CO 7 = 4 — > 3 ratio of 0. 15 is shown as a red line with circles. 



cupying the same volume, hence the line ratios are independent 
of beam dilution. For optically thin clumps (or for clumps with 
moderate opacity), the observed-to-model integrated intensity 
ratio gives the ratio of the beam-averaged column density to 
the clump intrinsic column density, assuming that the total in- 



tegrated intensity is proportional to the number of clumps, i.e. 
no velocity crowding is present. We therefore determined the 
beam filling by comparing the observed absolute intensity of a 
line with a constrained model. 



5 



Pineda, Mizuno, Rollig, Stutzki et al.: Submillimeter line emission from LMC 30 Dor 

19 



I 
O 

o 



o 




log [n(H 2 )/cm" J ] 

Fig. 4. Radiative transfer predictions and the observed absolute 
intensity of 13 CO J = 4 -> 3 for T kin = 160 K. The three red- 
dashed contours show the observed intensity and its value scaled 
up by factors of 8 and 16, respectively. The blue-dashed lines 
show the range of column and volume densities constrained by 
the analysis of the CO line ratios. The absolute intensity of 13 CO 
J = 4 — > 3 needs to be scaled up by a factor of about 16 to 
reproduce the constrained column and volume densities. 

Table 3. Parameters for different relevant molecular and atomic 
transitions. The listed critical densities were calculated for colli- 
sions with p-H2 and a kinetic temperature of 100K. The critical 
densities were calculated using the Einstein A- and collisional 
rate coefficients tak en from the Leiden Atomic and Molecular 
Database (LAMDA; ls"choier et alj|2005l) . 



Line 


Frequency/GHz 




%/cm 3 


12 CO J = 1 -> 


115.2712 


5.5 


2xl0 3 


12 CO J = 3 -» 2 


345.7959 


33.2 


4xl0 4 


12 CO J = 4 -> 3 


461.0408 


55.3 


9xl0 4 


12 CO J = 7 -» 6 


806.6517 


154.9 


4xl0 5 


12 CO J = 9 -> 8 


1036.912 


248.9 


lxlO 6 


13 CO J = 1 -» 


110.2013 


5.3 


2xl0 3 


13 CO J = 4 -» 3 


440.7654 


52.9 


8xl0 4 


[Cl] 3 P! ^ 3 P 


492.1607 


23.6 


lxlO 3 


[Ci] 3 P 2 -» 3 Pj 


809.3446 


62.5 


lxlO 3 


[Ol] 3 P, -+ 3 P 2 


4744.777 


227.7 


6xl0 5 


[Oi] 3 P -> 3 P, 


2060.069 


326.6 


7xl0 4 


[Cn] 2 P 3/2 ^ 2 P 1/2 


1900.5377 


91.2 


5xl0 3 



4.7. [CI]-emission 

In Table [2] we see that the ratio between the two [C 1] lines is 
0.92. We compare this ratio with predictions from the radia- 
tive transfer calculations in Figure |2] We present model line ra- 
tios for three different kinetic temperatures 80, 120, and 160 K, 
as a function of the [Ci] column density per velocity inter- 
val, Af[ci]/Av, and H2 volume density, n(H2). The coverage in 
Afcrj/Avis 10 15 -10 I9 cm- 2 andin«(H 2 )is 10 2 -10 6 cm" 3 . The 
observed ratio can only be reproduced by kinetic temperatures 
above 60 K and Hb densities of about 1 - 3 x 10 3 cm -3 . This Ho 




log [n(H 2 )/cm"' 



Fig. 5. Comparison between the observations and the predictions 
from the radiative transfer model for the 12 CO J = 1 — » and 
J = 4 — > 3 line ratios (black-solid contour) and the 13 CO to 
12 CO / = 1 — > line ratio (gray-dashed contour) for Tkin = 
160 K. The observed 13 CO to 12 CO / = 1 -> ratio of 0.08 is 
shown as a blue line with boxes while the 12 CO J = 1 — > to 
./ = 4 — > 3 ratio of 0.22 is shown as a red line with circles. 



density is similar to, or somewhat higher than, the critical den- 
sity of the [C 1] 3 Pi — > 3 Po line (~ 1 x 10 3 cm -3 ), and therefore its 
emission is close to being thermalized. From the analysis of the 
CO emission (see below) we found that the beam-filling factor 
in 30Dor-10 is about 1/16, so that the observed [Ci] 3 Pi - 3 P 
brightness of 0.58 K should be scaled up to 9.3 K. This bright- 
ness temperature corresponds to a [C 1] column density per ve- 
locity interval of about 1 x 10 17 cirr 2 /km s _1 . The multiplication 
by a filling factor is appropriate because, in the range of kinetic 
temperatures and H2 densities inferred above, the [C 1] 3 Pi - 3 Po 
line is likely optically thin. 



4.2. CO-emission 

We compare the observed 12 CO / = 7 -> 6/ 12 CO J = 4 -> 3 and 
13 CO 7 = 4^ 3/ 12 CO 7 = 4^3 ratios with predictions of the 
radiative transfer code in Figure [3] We assumed t he 12 CO/ 13 CO 
abund ance ratio of 35 derived in 30Dor-10 by Hei kkila et al.l 
dl999l) . We required that the volume density of the CO-emitting 
region is equal to, or somewhat higher than, that constrained for 
the [C 1] emission, because we expect that both [C 1] and mid-J 
CO lines are emitted from adjacent layers in a PDR-like struc- 
tural where the neutral carbon is located in an equal or lower- 
density region compared to the one at which CO is located, and 
that the cloud's volume density profile is smooth. In Figure [3] 
this condition is only valid for kinetic temperatures higher than 
120K for a volume density range between 3 x 10 3 - 10 4 cm~ 3 . 
The intersection of the two CO line ratios constrains the CO col- 
umn density per velocity interval to about 3 x 10 17 cirT 2 /km s~' 
at Tan =a 160 K. The derived kinetic temperature and H2 den- 



1 For an example of a typical temperature distribution and abundance 
st ructure of the main carbon sp ecies in PDRs, see e.g. Figures 7 and 9 
in lTielens & Hollenbachll 19851 . 
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sity are consistent with previ ously published excitat i on analyzes 
based on CO observations (|Johansson et al.l Il998t iKiml l2006t 
iMinamidani et alj|2008ll201 lb . 

Figure |4] shows that the model-predicted clump-averaged 
absolute line intensity of the 13 CO / = 4 — > 3 line, at the column 
and volume density regime constrained by the observed line ra- 
tios, is about 16 times higher than that observed at T^ n = 160 K. 
This implies a beam filling factor for CO of about 1/16, which 
corresponds to a typical clump size of 10". Recent ATCA obser- 
vations of HCO + revealed several compact clumps in 30Dor-10 
(Anderson et al. in preparation). The clumps have a typical an- 
gular size of 10", and are marginally resolved at the resolution 
of their observations, 6". More observations with higher angular 
resolution are required to confirm the observed clump sizes in 
30Dor-10. 

In Figure [5] we consider line ratios involving the 12 CO and 
13 CO / = 1 -> transitions for T kin =160 K. The 12 CO 
J - 1 — > to / = 4 — > 3 ratio of 0.21, however, sug- 
gests lower column and higher volume densities than those con- 
strained above. The low 12 CO J - I — > to ./ = 4 — > 3 ratio is 
due to the observed 12 CO J = 1 —> line that is too weak for 
the physical conditions we derived. The observed 13 CO/ 12 CO 
ratio is therefore also affected by this weaker 12 CO / = 1 — > 
emission. This may indicate an additional colder gas component. 
We were only able to reproduce absolute intensities of the 12 CO 
and 13 CO J = 1 -> lines using a model with the conditions 
obtained for the column and volume density derived for [Ci] 
(Section |4. \\ and kinetic temperatures of 7\i n =10-30K. This is 
consistent with observations of higher density tracers, suggest- 
ing a colde r (~50K) and denser (~ 10 s cm -3 ) gas component in 
30Dor-10 dHeikkilaet al.lfl999b . 

4.3. [Ci] and CO column densities 

Combining our estimate of the CO and C i column density per 
velocity interval, the observed line width of 9 km s _1 , and the 
beam filling factor of 1/16, we derive a beam-averaged column 
density of 1.7xl0 17 crrT 2 and 5.6xl0 16 crrT 2 for CO and Ci, 
respectively. Thus, the ratio of the column density of [Ci] to 
that of CO is ~3. Since the 12 CO and 13 CO J = 1 -> lines 
indicate an additional, cooler component, this ratio is an upper 
limit. Assuming that all gas-phase carbon is in the form of C i 
and CO in their respective line-emitting regions, we used a car- 
bon abundance relative to H 2 in the LMC of 1 .6 x 10~ 4 (de- 
rived from l"Cl/rHl=8x 1 5 : iDufour etaill 1 982b to obtain an H 2 
mass of 1625 M associated with CO and of 531 M Q associated 
with C 1. Thus, the total mass within the beam of our observa- 
tions is 2156 M . Note that our beam size of 38"FWHM cor- 
responds to 9.5 pc at the distance of the LMC and hence covers 
1/2 of the 30Dor-10 cloud's extent of about 20 pc. If the H 2 
column density is uniform over 30 Dor-10, the total mass of the 
cloud should be a factor of 4 larger than the mass derived within 
the beam of our observations, reflecting the difference in area. 
Thus, the total H 2 mass in 30 Dor-10 should be about 8.6xl0 4 
M Q . This derived mass agree s well with the CO-lum inosity mass 
of 8.5 x 10 4 M Q derived by Uohansson etaf] d!998l) assuming a 
CO-to-H 2 conversion factor of 4xl0 20 cm' 2 (K km s -1 ) -1 (e.g. 
IIsraelfl997l iHughes et al.ll20iot iPineda et al.ll2009l) . Note that 
some H 2 might be not traced by CO or [C 1] but by [C 11] (the 
"hidden H 2 gas"). Thi s can account for about 30-50% of the 
H 2 in the Galaxy (e.g iGrenier et al. 2005; Langer e t al.l 1201 Ct 
IVelusamv et al J 120101 : iPlanck Collaboration et al.l 1201 lb . and is 
expect ed to have an even larger mass in low-met allicity environ- 
ments dMadden et al.|[l997t IWoIfire et al.ll20"Iob . Therefore, the 



Table 4. Comparison between observed line integrated intensi- 
ties and the predictions from the KOSMA-t PDR model and 
the RADEX radiative transfer code. All intensities are in units 
of erg cm 2 s _1 sr _1 . The observed and RADEX intensities are 
corrected for beam dilution. The PDR model absolute intensities 
are calculated assuming a common FWHM source size of 90". 



Line 


Observed 


KOSMA 


-r" 


RADEX" 


12 CO J = 1 -> 


2.8xl0- 7 (0.2) 


8.7x10 


-7 


8.IXIO- 7 


12 CO J = 3 -> 2 


3.0xlO- 5 (0.3) 


3.4x10 


-5 


3.2xl0~ 5 


12 CO J = 4 -> 3 


7.8xl0- 5 (1.6) 


7.8x10 


-5 


6.8xl0~ 5 


12 CO J = 7 -» 6 


2.2xl0~ 4 (0.4) 


1.9x10 


-4 


2.1xl0" 4 


12 CO J = 9 -> 8 


<3.3xl0- 4 


9.4x10 


-5 


1.9xl0~ 4 


l3 CO J = 1 -> 


1.9xl0- 8 (0.3) 


6.8x10 


-8 


3.1xlO~ 8 


13 CO J = 4 -> 3 


1.0xlO- 5 (0.3) 


1.1x10 


-5 


1.0xl0~ 5 


[Ci] 3 Pi -> 3 P 


l.lxlO- 5 (0.3) 


1.0x10 


-5 


5.5xl0~ 6 


[Ci] 3 P 2 -^Pj 


4.1xl0- 5 (1.7) 


4.4x10 


-5 


3.5xl0 5 


[Cn] 2 P 3/2 ^ 2 P 1/2 


1.2xl0- 2 (0.01) 


2.4x10 


-2 


4.5xl0- 3 


(a) Best-fit model: 
Xo = 3160. 


Hens = 3.3 X 10 5 Cm 


~\ M em = 


1.2 x 10" Mo, 



(b) Best-fit model: n Sl = 10 4 cm" 3 , T m = 160 K, beam filling 1/16, 
N co = 1.7 x 10 17 cm^,iV c o = 5.6 x 10 16 cm 2 . 

mass estimate derived here should be considered a lower limit to 
the total H 2 mass. 

5. PDR-model analysis 

In the following we investigate whether our observations toward 
30 Dor- 10 can be explained in terms of an ensemble of PDRs dis- 
tributed in a cl umpy interstellar medium. We used the KOSMA- 
t PDR model dStorzer et alHl996l: iRollig et alj|2006b to model 
an ensemble of spherical clumps with a power-law mass spec- 
trum dN/dM cc M ~ L8 and a mass-size relation M oc r 23 
(ICubick et al.ll2008b . The model provides a self-consistent so- 
lution of the chemistry and thermal balance. The free parame- 
ters of the clumpy PDR ensemble are (1) average ensemble total 
H volume density («H+2nH 2 X "ens, (2) ensemble mass, Mens, (3) 
strength of the FU V radiation field in units of the iDrainel (11978b 
field, xo, and (4) the minimum and maximum mass of the clump 
ensemble, [M llun ,M max ]. We fix the metallicity in the model to 
Z = 0.4. We attempted to fit the observed lines with a model 
with Z Q = 1, but this solar metallicity model could not repro- 
duce the mid-/ CO observations with a^- 2 about factor of 10 
larger than for the low-metallicity model. Note that we assumed 
that the clumps do not spatially overlap, therefore we considered 
optical depth effects only within individual clumps. We also as- 
sumed that the FUV field is uniform within the beam of the ob- 
servations. 

We constrained the model with the absolute intensities of 
all five lines observed with NANTEN2 together with the 12 CO 
J — 3 — > 2 and [Cn] 158/im lines taken from the literature 
(Section 13.2b . Fitting absolute line intensities instead of line 
ratios allowed us to study the internal structure of the cloud 
(dumpiness) via the area and volume filling factors, and it adds 
the cloud mass and size as model parameters. The source size 
was fixed to be 90" (Section [3.1t . We did not include the 12 CO 
and 13 CO J — 1 — » transitions in the model fit because these 
transitions are likely affected by optical depth effects between 
clumps that are not accounted for by the model. We still com- 
pared the intensities of the 12 CO and I3 CO / = 1 — > transi- 
tions with the predictions of the PDR model and discuss below 
the possible reason of the discrepancy between model and obser- 



7 



Pineda, Mizuno, Rollig, Stutzki et al.: Submillimeter line emission from LMC 30 Dor 



vations. We used simulated annealing to find the optimum model 
parameter combination that predicts absolute line intensities that 
match those observed. 




4 4.5 5 5.5 6 

log [n ens /cm" 3 ] 



Fig. 6. Comparison between the different observed lines with the 
predictions from the clumpy-PDR model for a given H2 volume 
density and FUV radiation field. In all panels, the solid lines rep- 
resent the observed values while the dashed and dotted lines rep- 
resent their lower and upper error bars, respectively. The gray- 
scale represents the value of x 2 relative to its minimum value, 
X 2 - ■ 

A mm 

The best-fit model results suggest an average ensemble to- 
tal H volume density n eBS = 3.3 x 10 5 cm~ 3 , an ensemble mass 



M ens = 1 .2x 10 5 Mq. The minimum and maximum masses of the 
clump ensemble are M min = 1 x 1(T 3 M Q and M max = 0.1 M Q . 
Note that « ens corresponds to the average total H volume density 
(i.e. «h+2«Ho)- Ignoring the contribution from atomic H, the av- 
erage ensemble H2 volume density is about 6 x 10 4 cirT 3 , consis- 
tent with that derived with the excitation analysis in Section [4] 
The ensemble mass is also consistent with that derived for the 
whole extent of the 30Dor-10 cloud in Section 1431 A compar- 
ison between the observations, the predictions from the PDR 
model, and the results of the excitation analysis (Section |4]i is 
listed in Table g] 

The area filling factor of the clump ensemble is defined as the 
sum of all clump projected areas over the cloud solid angle (as- 
sumed to be 90"). Similarly, the volume filling factor is defined 
as the sum of all clump volumes over the volume of a sphere with 
90" projected diameter. The model fit gives an area filling fac- 
tor of 2.94 and, assuming a distance of 50kpc, a volume filling 
factor of 0.001. These area and volume filling factors imply that 
the model requires many low-mass clump s to reproduce the ob - 
served line intensities. As discussed by e.g. lKramer et al.ld2008l) . 
these low-mass clumps would require very high external ther- 
mal pressures to be confined and might be considered as tran- 
sient features of the turbulent gas and should evaporate on short 
timescales. Note that the area filling factor defined here is not 
comparable with the beam filling factor derived in Section 14.21 
because the PDR model does not predict whether these small 
clumps are clustered or dispersed within the cloud. 

In Figure|6]we show the dependence of the different observed 
transitions on different H2 volume densities and FUV radiation 
fields. Solutions in the range « ens = 1 X 10 5 - 1 x 10 5 6 crrT 3 and 
Xo = 400 - 6300 are within three times the minimum;^ 2 • = 59. 
Note that this high value ofx 2 is dominated by the inability of the 
PDR model to reproduce the [C 11] emission. If we exclude the 
[Cn] line in the calculation of^ 2 , we obtain a much lower x 2 lin 
of 0.01, without significant changes in the constrained physical 
conditions. 

To illustrate the importance of the observed line intensities 
in determining the value of xo and « ens , we show in Figure [7] 
the predicted model line intensities as a function of n ens for a 
fixed value of;^o=3980, and as a function of^o for a fixed value 
of « ens =3.9xl0 5 crrT 3 . The fixed values of^o and « ens were se- 
lected to be the closest point in the model grid to the solution 
found using the simulated annealing technique. The intensity of 
the mid-/ 12 CO and 13 CO are very sensitive to « ens and drive 
the value of this parameter. The [C 1] and mid-/ CO lines are 
sensitive to the FUV radiation field but to lesser extent than to 
«ens- Consequently, as also seen in Figure[6] the fit to the average 
ensemble volume density is better than that for the FUV radia- 
tion field. The decreasing value of the [C 1] and mid-/ CO lines 
and the increase of value of the [C 11] intensity with xo is related 
to the shift of the C + /C°/CO transition deeper into the cloud as 
Xo increases, reducing the column densit y of CO and increas - 
ing column density of [Cn]. As shown by ICubick et al.l d2008l) . 
the emission from CO transitions with / > 5 is mostly arising 
from low-mass clumps, while lower-/ CO and [C 1] transitions 
arise from clumps with larger masses. Therefore, the observed 
I2 CO / = 7 — > 6 transition is very sensitive to the lower mass 
limit and drives the value of this model parameter. For the upper 
mass limit, the dependence on the observed lines is weaker, and 
varying the upper mass limit to higher values does not have a 
significant effect on the model fit. 

In Figure [6] and Table |4] (see also Figure [8), we can see that 
the best-fit model agrees well with observations of mid-/ CO 
and [Ci] but it overestimates the [Cn], 12 CO / = 1 -» 0, and 
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3 P 
[CI] ^P^fPn 




log [n ens /crrf 



Fig. 7. Model-predicted integrated intensities the spectral lines 
used to constrain the PDR model as a function of xo for a fixed 
value of rc ens {upper panel) and of n ens for a fixed value of xo 
{lower panel). 



13 CO 7 = 1 — > line intensities. To reproduce the lines ob- 
served by NANTEN2, the model requires clumps with volume 
densities that are much higher than the critical density of the 
[Cn] line (5 x 10 3 cm' 3 ; see Table |3j, resulting in a predicted 
[Cn] intensity that is higher than observed. A model with a 
lower volume density might reproduce the [C n] line better, but 
it would underestimate the mid-7 and [Ci] transitions. Another 
complication is that, as shown in Figure [T] we can expect that 
the FUV field is not uniform within 30Dor-10, because parts 
of the cloud are closer to R136 than others, while the model 
assumes a uniform field. Note that we are only considering a 
single-component model in this analysis. It might be possible 
to find a better fit to the lines if more gas components were 
added, but more data (e.g. [Oi] lines) would be needed to ob- 
tain a good constraint to the model. Note that, as mentioned 



above, the PDR model neglects any mutual absorption/shielding 
between clumps. Since the 12 CO 7 — 1 — > is optically thick 
and I3 CO 7=1 — > is likely affected by optical depth effects, 
radiation from background clumps should be (partly) shielded 
by foreground clumps. Because the model ignores this effect, it 
overestimates the intensity of these transitions. This effect might 
also partly be the reason for the overestimation of the [C n] line, 
because this line might be also affected by optical depth effects 
between clumps. 

Our determination of the strength of the FUV field agrees 
excellently wi t h the independent estimation presented by 
iPoglitsch et ail d 19951) of xo = 3100-3600. The constraints on 
the H2 volume density and FUV radiation field for 30 Dor are 
also consi stent with previous determinations b a sed on PDR 
modeling (iPak et all U998T: iKaufman et ail 1 19991 : iBolatto et alj 
Il999t lR611ig et al. 2006). A crucial test to the suggested clumpy 
structure of the LMC's ISM will be high angular resolution 
imaging of low- and mid-7 CO and [Ci] transitions with 
ALMA. 



6. Discussion 

6.1. Comparison with LMC-N159W 

The 30Dor-10 region is located at a distance of ~ 20 pc from 
the center of the R136 star cluster, hence this region is exposed 
to an extreme FUV radiation field. It is of interest to compare its 
physical properties with other regions in the LMC that have more 
moderate ambient radiation fields. Since there ar e no significan t 
spatial variations of the metallicity in the LMC (lDufourlll984l) . 
this comparison will help us to understand the effect of the FUV 
radiation field on low-metallicity environments. 

In iPineda et all d2008l) we made a similar analysis of the 
same set of spectral lines as analyzed here in the LMC-N159W 
region. We constrained the H2 volume density to be ~ 10 4 cm -3 , a 
kinetic temperature of ~80 K and equal CO and C 1 column den- 
sities of ~1.6xl0 17 cirT 2 . The beam filling in N159W is about 
1/6. Modeling of the photon-dominated region suggests that the 
strength of the FUV field in N159W is about ^ =220, while we 
find^o = 3160 in 30Dor-10 (Section^. 

The 30Dor-10 region is warmer than N159W, with a ki- 
netic temperature of about 160 K and 80 K, respectively. This 
higher temperature is indicated by the higher 12 CO 7 = 7 — » 6 
to 7 = 4 — » 3 ratio in 30Dor-10. This warming can be un- 
derstood because the gas heating is likely dominated by photo- 
electric heating, which d epends on the strength of the FUV field 
dBakes & Tielensl 1 1 998b . The higher temperatures in 30Dor- 
10 do not result in higher brightness temperatures compared 
with N159W because of the reduced beam filling (1/16). The 
column/volume densities, however, are similar or somewhat 
higher in 30Dor-10 compared with N159W, suggesting that CO 
clumps in the former region are smaller. These smaller clump s 
can be interpreted as the effect of the enhanced CO photo- 
destruction due to the strong FUV radiation field. This enhanced 
photo-destruction does not result in a higher abundance of C 1 
relative to CO, becuase the N{CO)/N{C 1) column density ratio is 
even higher in 30 Dor-10, implying that the FUV field moves the 
C + /C/CO transition layer to locations deeper into the clumps, but 
does not increase the thickness of this layer. Accordingly, most 
of the carbon in 30 Dor-10 should be in the f orm of C + . This ef - 
fect has been predicted by PDR models (e.g. lRollig et alj|2006l) . 
Note that we cannot eliminate other mechanisms that might de- 
stroy the molecular cloud, such as stellar winds from massive 
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stars. However, the N159W should also be influenced by these 
effects. 

We determined the column density of C + using the 
integrated intensity observed by iBoreiko & Betzl (Il99ll) of 
7.7xlO~ 4 ergcirT 2 s sr _1 . This quantity is related to the C + col- 
umn density, assuming a kinetic temperature larger than 91 K 
and a volume density higher than « cr = 5 x 10 3 c m -3 , by N(C + )= 
6.4 x 10 20 /([Cii]) cm-^ergcnrVsr- 1 )- 1 dCrawford et al.l 
fl985l) . With that we obtain a column density of C + of 
4.96xl0 17 cirT 2 and therefore the total carbon abundance in 
30Dor-10 is distributed as JV C + : N c : N C0 =69%:8%:23%. 
For comparison, the distribution of carbon in the N159W re- 
gion is 46%:27% :27% while in the Galactic PDR DR21(OH) it 
is 3%:10%:87% dJakob et alJIiOOl . As in N159W, the 30Dor- 
10 region has an increased fraction of the gas-phase carbon in 
the form of C + compared with DR21(OH). 

6.2. CO spectral-line energy distribution 

We combined our observed mid-7 12 CO observations with pre- 
viously observed transitions to derive the CO spectral-line en- 
ergy distribution (SLED) in 30Dor-10 and N159W. We present 
the CO SLEDs in Figure [8] together with the predictions from 
our excitation analyses and PDR modeling i n both regions. We 
includ e observations of 12 CO / = 3 — > 2 from lMinamidani et al.l 
(|200 8j) and upper limits to t he intensity of the / = 9 — > 8 line 
from lBoreiko & Betzl d 199 lb . We corrected these line intensities, 
together with that of the J = 7 — > 6 transition, to match the 
intensity corresponding to a 38" beam (see Section[3j. Single- 
component models agree well with the observed CO transitions. 
But the model for 30Dor-10 still predicts a line intensity of the 
/ = 1 — > line that is about a factor of 3 higher than ob- 
served. The 12 CO SLED peaks at / = 8 -> 7 for 30Dor-10 
and J = 6 -» 5 for N159W. The different CO excitation in both 
regions reflects the difference in their kinetic temperature, which 
in turn is due to the difference in the FUV radiation field between 
both regions. 

The CO SLED predicted by the clumpy PDR model for 
30Dor-10 is similar to that predicted by the excitation analysis 
for low-/ transitions, but differs significantly for high-/ tran- 
sitions, with the PDR model predicting much higher intensities 
for / > 6. For / > 14, however, the PDR model predicts almost 
constant intensities, while the intensity predicted by the excita- 
tion model drops. The clumpy PDR model predicts a peak of 
the CO SLED at a lower 12 CO transition (/ = 6 — » 5) compared 
with the prediction from the excitation analysis (/ = 8 — » 7). For 
N159W the clumpy PDR model CO SLED also predicts lower 
intensities for high-/ transitions but the difference is smaller 
than for 30Dor-10. The peak of the CO SLED still coincides 
with that predicted by the excitation analysis ( 12 CO / = 6 — > 5). 
More observations of high-/ CO transitions with e.g. Herschel 
or SOFIA in 30 Dor-10 and N159W are required to test the pre- 
dictions of the PDR modeling and excitation analysis for high-/ 
CO transitions. 

In Figure|8]we also show the CO SLED for two positions in 
the Carina region, Carina-N and Carina-S, whose physical con- 
dit ions have been const rained with the KOSMA-t PDR model 
by |Krameretai](l2008l) using the same set of line as presented 
here. The constrained values of (xo, « ens ) are (3 100, 2x 10 5 cm -3 ) 
for Carina-N and (310, 2 x 10 5 cm 4 ) for Carina-S. Therefore, 
both Carina-N and Carina-S have similar physical conditions 
compared with 30 Dor-10 and N159W, respectively, but solar 
metallicity. The metallicity, however, seems not to affect the 
shape of the CO SLED because for both pairs of sources the CO 



spectral energy distribution is consistent within the small differ- 
ences in the physical conditions and uncertainties of the models. 

The total CO luminosity (Leo) predicted by the radia- 
tive transfer code in an area defined by the 38" beam of our 
NANTEN2 observations (corresponding to 9.5 pc at a distance 
of 50kpc) is l.lxlO 36 erg s" 1 for 30 Dor-10 and 4.3xl0 35 erg 
s _1 for N159W. The clumpy PDR model predicts a somewhat 
lower Leo of 8.2xl0 35 erg s 1 for 30 Dor-10 and 2.5xl0 35 erg 
s -1 for N159W. From the observations, we also estimated the 
total [C n] and [C i] luminosity (ignoring the much weaker [C i] 
3 P 2 -> 3 Po transition) of 9.9xl0 36 and 4.2xl0 34 erg s _1 , respec- 
tively for 30 Dor-10 and 2.6xl0 36 and 4.9xl0 34 erg s" 1 , respec- 
tively, for N159W. The relative fraction of the luminosity of 
these species is L [C ii] : L [C rj : L CO =90%:0.4%:9.6% for 30 Dor- 
10 and 84.4%:1.6%:14 for N159W, when using L co predicted 
by the radiative transfer code, and 92%:0.4%:7.6% for 30 Dor- 
10 and 90%:1.7%:8.3% for N159W, when using L co predicted 
by the clumpy PDR model. We also estimates a total [Oi] lumi- 
nosity of 7.4xl0 36 erg s" 1 from observati ons of the [Oi] 63yum 
and 146/^m lines bv lPoglitsch et al.l d 19951) in 30 Dor. These esti- 
mated luminosities suggest that [C n] and [O i] contribute equally 
to the gas cooling in 30Dor-10. This significant contribution of 
[Oi] to the line cooling confirms the warm and dense gas con- 
ditions derived with the excitation analysis and PDR modeling. 
Note that we cannot rule out the possibility of an additional com- 
ponent powering even higher CO transitions (/ > 10) in 30 Dor- 
10 and N159W. Therefore, our estimate of the total CO luminos- 
ity might be a lower limit. As mentioned above more observa- 
tions will allow us to determine the contribution of high-/ CO 
transitions to the total CO luminosity. 

7. Summary and conclusions 

We derived the physical conditions of the line-emitting gas in 
the 30 Dor-10 region in the LMC. We compared an excita- 
tion/radiative transfer code and a PDR model with NANTEN2 
observations of the 12 CO / = 4 -> 3, / = 7 -> 6, and 13 CO 
/ = 4 -> 3 rotational and [Ci] 3 Pi- 3 P and 3 P 2 - 3 Pi fine- 
structure transitions. Our results can be summarized as follows: 

- The analysis of the excitation conditions for both the CO 
submm-lines and the [C i] fine structure lines shows temper- 
atures of about T = 160 K and densities of about «h, = 10 4 
cm -3 for the emitting gas. 

- We found that 30 Dor-10 is warmer and has a lower beam 
filling factor compared to N159W. This difference might be 
the result of the FUV radiation field heating the gas and pho- 
todissociating CO molecules. This, however, does not result 
in an enhanced C abundance relative to CO: we obtained a 
N(CO)/N(Ci) ratio of 3 in 30 Dor-10, which is higher than 
that found in N159W of iV(CO) /N(C i) equal to unity. 

- We derived the CO spectral-line energy distribution of 
30 Dor-10 and N159W. Considering the excitation condi- 
tions constrained using line ratios, we found that this dis- 
tribution peaks at 12 CO / = 8 -> 7 for 30 Dor-10 and at 
/ = 6 — » 5 in N159W. The clumpy PDR model, however, 
predicts the peak of the CO spectral-line energy distribution 
at / = 6 -> 5 for both 30 Dor-10 and N159W. 

- We compared our observations with the results of a clumpy 
PDR model. The model that best reproduces the observed 
absolute integrated intensities has an average ensemble H 
density of n ens ~ 10 5 cirT 3 , a total mass of the ensem- 
ble, M ens ~ 10 5 M Q , and a strength of the FUV field of 
Xo ~ 3100. The constraints on the PL volume density and 
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Fig. 8. Spectral-line energy distribution (SLED) for several millimeter and submillimeter lines toward 30Dor-10 (left), N159W 
(center), and Carina-N and -S (right). In the left and central panel, the solid lines represent the predic tions for the 12 CO and 
13 CO SLED from the excitation analysis presented in Section |]for 30Dor-10 and in lPineda et al.l(l2008l) for N159W. The gray- 
dashed lines represent the 12 CO and 13 CO SLED predicted by the clumpy PDR model fits. The excitation analysis and PDR model 
predictions for the [C i] and [C n] lines are represented by triangles and boxes, respectively. In the right panel we show the 12 CO 
SLED for the Carina-N and Carina-S regions presented by iKramer et al.l d2008l) using the same clumpy PDR model as used in this 
paper. The physical conditions derived from the PDR model are similar for 30Dor-10 and Carina-N as well as for N159W and 
Carina-S, but they have different metallicity (Carina Z=Z ; LMC Z=0.4 Z ). 



the strength FUV radiation field agree well with independent 
determinations and previous PDR modeling results. 

The 30Dorarus region is the closest example of vigorous 
star-formation taking place in a low-metallicity environment, 
which is thought to be common in the early universe. Our re- 
sults suggest that the star-forming ISM in 30Dor-10, and per- 
haps in starburst galaxies in the early universe, is dense, and 
warm, and clumpy. The clumpy structure of the 30 Dor- 10 gas 
is a result of the extreme FUV environment to which the low- 
metallicity gas is exposed. Our results also highlight the impor- 
tance of ionized carbon in tracing the total gas mass and, to- 
gether with atomic oxygen, in regulating the thermal balance of 
the gas. Future observations of fine-structure transitions of [C n], 
[O i], [N n], and [C i] and rotational transitions of CO in the local 
Universe are essential for our understanding of star-formation 
in low-metallicity environments and are an important tool for 
the interpretation of observations of high-redshift galaxies with 
ALMA. 
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